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CARLEMAN ESTIMATE FOR BIOT CONSOLIDATION SYSTEM IN PORO-ELASTICITY 
AND APPLICATION TO INVERSE PROBLEMS 

M. BELLASSOUED AND B. RIAHI 


Abstract. In this paper, we consider a coupled system of mixed hyperbolic-parabolic type which describes 
the Biot consolidation model in poro-elasticity. We establish a local Carleman estimate for Biot consilidation 
system. Using this estimate, we prove the uniqueness and a Holder stability in determining on the one hand a 
physical parameter arising in connection with secondary consolidation effects A* and on the other hand the two 
spatially varying density by a single measurement of solution over cu x (0, T), where T > 0 is a sufficiently 
large time and a suitable subbdomain cj satisfying dcj D dfl. 


1. Introduction 

Let us consider an open and bounded domain O of with boundary F = dfi. Given T > 0, the 
Biot consolidation model in poro-elasticity in which we are interested is the following 

I Uit - - V(A*(x)divut) + pi(x)V 6 '= f in Q = 0 x (0,T), 

Ot - A0 + g 2 (x)divut = g in Q, 

with Dirichlet boundary condition 

(1.2) u(x,t)=0, d(x,t) = 0, onS = rx(0,T), 

where the .t stands for the time derivative, V = (cli, 92 , ^ 3 ), and a is the elliptic second order linear 
differential operator given by 

A^^av(x) = gAv{x) + ig + A)(V(divv(x))) 

(1.3) + divv(x)VA(x) + (Vv + (Vv)^)V/r(x), x G 12, 

for V = {vi,V 2 ,vs)'^, where denotes the transpose of matrices. Throughout this paper, t and x = 
(xi, X 2 , X 3 ) denote the time variable and the spatial variable respectively, and u = (ui, tt 2 , 1 x 3 )^ denotes the 
displacement at the location x and the time t, and 9 = 9{x, t), the temperature, is a scalar function, 5 is a 
heat source. We will assume that the Lame parameters g, X G C^(12), satisfy 

/r(x) > 0, A(x) + 2|u(x) >0, Vx G 12. 

We can prove (e.g., |[T] |371) that the system (11.11) possesses a unique solution (u, 9) with suitable initial 
values. 

The main subject of this paper is the inverse problem of determining, not only, a physical parameter 
arising in connection with secondary consolidation effects A* but also, the two spatially varying density 
Pi and Q 2 , in the Biot consolidation model in poro-elasticity, uniquely from observed data of displacement 
vector u and the temperature 0 on a suitable subdomain a; C 12 and the observation data of u and 9 at given 
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a suitable time Iq. Such kinds of observation data are similar to those considered in (e.g. |6I), which are 

typical for obtaining the corresponding stability results by the Carleman estimates. 

I. 1. Inverse Problem. Let cj C be a given arbitrarily subdomain such that Stu D F. The condition 
duj D r, is used to deal with the lack of a divu boundary condition. More precisely, to overcome the 
difficulty due to the strong coupling in the system (11.11) . we have to find the Carleman estimate for divergence 
of the first equation of the system (11.11) . But we do not know the boundary data of divu. Due to this reason, 
we convert the estimates of divu on the boundary F to the estimates in w, a neighborhood of the whole 
boundary. 

The goal of this paper is to derive a Holder stability and the uniqueness, by measurements 

uLx(0,t), and 9{x,to) x G D. 

The key ingredient in our argument is an L^-weighted inequality of Carleman type for consolidation Biot’s 
system. This coupled mixed hyperbolic-parabolic system can describe the phenomena arising when a soil 
is submitted to a load as well as the ultrasonic propagation in fluid-saturated parous media like cancellous 
bone. The displacement vector field of the system, denoted by, u, satisfies the conservation of momentum 
while the fluid pressure, 9, satisfies a diffusion equation. 

K.Terzaghi 1381 was the first in interesting in the consolidation phenomena arising in porous media under 
a load. He showed the similarities between this phenomena and the exit of a flow out of a porous media, 
that contributed to model fluid flows in saturated deformable porous media as a coupled flow-deformation 
process. 

Later, M. A. Biot [H |9l [TOl [TT] [121 studied these problems assuming that the continuum mechanics 
laws are applicable. He develop thus the now classical theory of poro-elasticity and proved that the linear 
theory of consolidation could be established by using the Darcy law for laminar flows combined with the 
momentum balance equations with Hooke law for elastic deformations. 

In the particular case when the secondary consolidation term A* = 0 in system (11.11) . we are interested in 
the thermoelastic model: 

I putt{x,t) - A^^xu{x,t) + gV9{x,t) = ¥{x,t) in(5 = Dx[0,T] 

I co9t{x,t) — A9{x,t) + gdivut = h{x,t) in Q 

Bellassoued and Yamamoto Q established Carleman estimates with second large parameter for a coupled 
parabolic-hyperbolic system, a thermoelastic plate system and a thermoelasticity system with residual stress. 
According to the linear theory of thermoelasticity, Bellasoued and Yamamoto [4] consider a bounded and 
isotropic body whose mechanical behaviour is described by the Lame system coupled with the heat equation. 
Assuming the null surface displacement on the whole boundary. They prove a Holder stability estimate for 
the inverse problem of determining the heat source only by observation of surface traction on a suitable 
subdomain along a sufficiently large time interval using Carleman estimate for thermoelasticity system. 

Then, in a thermoelastic model, B.Wu and J.Liu [!6| study the inverse problem of determining two spa- 
ciahy varying coefficients with the following observation data: displacement in a subdomain ui satisfying 
duj D 5D along a sufficiently large time interval, both displacement and temperature at a suitable time over 
the whole spatial domain. Based on a Carleman estimate on the hyperbolic-parabolic system, B.Wu and 

J. Liu prove the lipschitz stability and the uniqueness for this inverse problem under some a priori informa¬ 
tion. 
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Our method is based on the tool of Carleman estimates, which was originally introduced in the field of co¬ 
efficient inverse problems for hyperbolic, parabolic and elliptic equations with the lateral data by Bukhgeim 
and Klibanov ifT^ . 

For the formulation of inverse problems with a finite number of observations, Bukhgeim and Klibanov 
|[T5]| proposed a remarkable method based on a Carleman estimate and established the uniqueness for inverse 
problems of determining spatially varying coefficients for scalar partial differential equations. See also 
Bellassoued and Yamamoto [31, L4], El Bukhgeim llT3l . Bukhgeim, Cheng, Isakov and Yamamoto iLTdll . 
Imanuvilov and Yamamoto ||23| . Isakov |[26l . Khaidarov ll30l . Klibanov OTI . 13^ . Klibanov and Timonov 
1341, Klibanov and Yamamoto l35l . 

However, to the best authors’s knowledge, the inverse problem for the Biot consolidation model in poro- 
elasticity, especially in the cases of multiple coefficients, have not been studied thoroughly yet. The main 
difficulties for these inverse problems come from: 

(i) In the Carleman estimates we should choose the same weight function to deal with equations of 
different types, namely parabolic equation and a strongly damped hyperbolic equation. 

(ii) The interaction between the two equations due to the strong coupling of displacement and tempera¬ 
ture requires much more complicated mathematical analysis. 

There are not many works concerning Carleman estimates for strongly coupled systems of partial differ¬ 
ential equations where the principal parts are coupled. 

1.2. Notations and statement of main results. In order to formulate our results, we need to introduce the 
following assumptions. 

Assumption (A.l): Let xq G M^\n. Then we introduce the conditions on the scalar functions ^ and 2/i-|-A: 

(1.5) /^) A G /u(x) > /To > 0) 2^(x) -|- A(x) > x £ Q 

and there exist ro G (0, /tq) and ri G (0, /ri) such that 

( 1 . 6 ) ^ |V(log//)| |x - xol < 1 - —, ^ |V(log(2^-f A))| |x - xol < 1 - — x £ Q. 

2 /io 2 /ii 

Assumption (A.2): We assume that the solution (u, 9) satisfies fhe a priori boundedeness: 

(1-7) l|u||vy7,oo(-Q) -t- ||^||vi/5,oo(q) < Mo, 

for some given positive consfanf Mq. 

Throughouf fhis paper, lef consider fhe admissible sef for fixed sufficienfly smoofh functions ai and a 2 on 

r 

( 1 . 8 ) A = {A* GC 2 (n),||A*||j 2 (j^) <m,X* = ai,VA* = aa onr,A*(x) > Aq > 0 in 0} 

for fixed m > 0. Moreover, for fixed sufficienfly smoofh functions oq and a on T and posifive consfanfs 
M, go and e, we sef 

(1.9) B = {(pi, Pa) G (C^(n))^, ||pi||^ 2 (Q) + ||£'2|Ic2(q) ^ 

Pi = ao, Vpi = a on r, pi(x) > po > 0 in cj}. 

For (v, y) £ x H^{Q), we denote fhe following norms af times to by: 

3 

Nto{^,y) = || 2 /(-,to)|li 2 (Q) +^||div5^v(-,to)|li2(n) 

j=0 


( 1 . 10 ) 
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and 
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(1-11) ^to(v,y) = lly(->io)ll^3(n) + ^ ||div54V(-,to)||5f3(n) 

i=i 

Finally, let to E (0, T) such that 

(1.12) -|x — XqI < mm{to, r — to}- 

x&n 

Throughout this paper, we assume that /r and 2/i + A satisfy (11.51) and (u, 9) satisfies the a priori boundede- 
ness (11.71) . 

Theorem 1.1. We assume that 

and |divut(x,to)| > e. 

Let A*, A* E A, to E (0, T) such that ( 17.721) and let assumptions {A.l) and {A.2) be held. Then, there exist 
constants C = C{Ll, T, to, Mo) > 0 and k E (0,1) such that 

(1.13) ||A* - A*||^.(^) < C (||u - n||^4(,,(o,r)) + ^t^iu -5,0-0))". 

By Theorem ll.il we can readily derive the uniqueness in the inverse problem 

Corollary 1.2. Under the same assumptions as in Theorem U . l\ and if 

u(x,t) = n{x,t), {x,t) E a; X (0,r), 

u(.,to) = u(.,to), d{-Ao) = 0{.,to) X gLI. 

Then, A* = A*, in LI. 

Theorem 1.3. We assume that 

|V0(x,to) • {x - xo)| > e. 

Let (pi, P 2 ), (qi, Q 2 ) G B, to G (0, T) such that ( 17.721) and let assumptions (^.1) and (^.2) be held. Then, 
there exist constants C = C{Ll, T, to, Mo) > 0 and 5 G (0,1) such that 

(1.14) ||pi - Pill^i(o) + II02 - 02|li2(f^) < C (^115 - u||^7(^x(0,t)) + ^io(u - wd- ^)) • 

By Theorem 1 1.31 we can readily derive the uniqueness in the inverse problem 
Corollary 1.4. Under the same assumptions as in Theorem \1.3\ and if 

u{x,t) = u(x,t), (x,t) G cc X (0,T), 

u(.,to) = u(.,to), ^(-Ao) = Of,to) X gLI. 

Then, Qi = Qi and Q 2 = Q 2 in Al. 

The remainder of the paper is organized as follows. In section 2, we give a Carleman estimate for Biot’s 
consolidation system. In section 3, we prove Theorem 11.11 and Theorem 11. 31 
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2. Carleman estimate for Biot’s system 

In this section we will prove Carleman estimates for the Biot’s consolidation system. In order to formulate 
our Carleman type estimates we introduce some notations. Let i? : —)■ M be the strictly convex function 
given by 

(2.1) 'd{x) = \x — xo\^ , X G n, 

where xq ^ fi. Let us define to by 

(2.2) min{fo, (T - to)^} > ( mc^i)(x) ) , 

V x£Q. j 

and let 

(2.3) ='d{x) - j3 {{t-tof - . 

Fix 5 > 0 and /3 > 0 such that 

(2.4) /3min{fg, (T — fo)^} > m^i)(x) +5, 0 < /3 < tq. 

Then the function il){x,t) verifies fhe following properties 

(2.5) 0)</3M- 6, 'ipix, T)</3M- 5, for all x G 
fhere exisfs e G (0, T/4) such fhaf 

(2.6) maxi/)(x, f) </3M —for all f G (0,2e) U (T — 2e, T), 

xen 2 

and 


(2.7) 


min'0(x,fo) > I3M. 


We define now fhe weighf function y? : 12 x M —)• M by 

(2.8) (/p(x,f) = 7>0, 

where 7 is a large paramefer selecfed in fhe following and lef 

(2.9) a = a{x,t) = s'yif{x,t). 
We use usual function spaces, H^{Q), and 

^ ^1(0, T; L^iSl)) n L2(0, T; 
Lef (v, y) a solufion of fhe linear Biol consolidalion system 


( 2 . 10 ) 


Vtt{x,t) - A^^Av(x,f) - V(A*divvt(x,f)) + giVy{x,t) = f(x,f) 
yt{x,t) - Ay{x,t) + P 2 divvt(x,f) = h{x,t) 


in (5 = 12 X (0, T) 
in Q 


such lhal 


(2 11) Supp(v(-,t)) C 12, Supp(y(-,t)) C 12, for all t G (0,T). 

c^v(x,0) = c2/v(x,T) = 0, y{x,0) = y{x,T) = 0 for all x G 12, j = 0,1. 

A Carleman eslimale is an inequality for a solufion lo a partial differential equation wilh weighted L^- 
norm and is effeclively applied for proving fhe unique continuation for a partial differential equation wilh 
non-analylic coefficienls. As a pioneering work concerning a Carleman eslimale, we can refer lo Carleman’s 
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paper |[T^ where what is called a Carleman estimate is proved and applied it for proving the uniqueness in 
the Cauchy problem for a two-dimensional elliptic equation. 

The following theorem is a weighted Carleman estimate with second large parameter for Biot’s consoli¬ 
dation system ( 12 . 101 ) with assumption ( 12 . 11 b . 

Theorem 2.1 (Carleman estimate for Biot’s consolidation system). There exist two constants 7 * > 0 and 
C > 0 such that for any 7 > 7 *, there exists s* = s*( 7 ) > 0 such that the following estimate holds 

(2.12) C J I |v|^-|-|divv|^-|-Idivv^l^-|-I Vdivv|^-|-(T iVdivv^p 

+ \yfdt < J (|fp-|-| Vf |^-|- 7 ~^ct |/i|^)e^*‘^dx dt 

-I- j (^7~^ |Adivv(x,fo)P + |divvt(x,fo)P + |divv(x,fo)P + |Vdivv(x, fo)|^) e^^'^dxdf 

for any solution (v, y) G H^(Q) x to problem ( I2.i0l ) which satisfy ( I2.iil ) and any s > s*. 

In order to prove Theorem l2. 11 we need a parabolic, strongly damped hyperbolic and hyperbolic Carleman 
estimates with second large parameter. 


2.1. Carleman estimate for parabolic equation. We consider the second order parabolic operator dt — A. 
As for Carleman estimate for parabolic equations with singular weight function, we can refer to Fursikov 
and Imanuvilov ifTTll . Imanuvilov lfT9l . Imanuvilov and Yamamoto 1201 . Here we give a Carleman estimate 
for parabolic equations with the regular weight function p. 

The following parabolic Carleman estimates holds 


Lemma 2.2. Let /c G N U {0}. There exist three positive constants 7 *, s* and C such that, for any 7 > 7 * 
and any s > s*, the following inequality holds: 

(2.13) Cj f + + 

\ a \=2 


< [ |(yt — div(A*Vy))|^ e^^'^dx fit 

Jq 


for any y G with compact support in Q. 

As for the proof, we can refer to Bellassoued and Yamamoto |j3l. 


2.2. Carleman estimate for strongly damped hyperbolic equation. In this subsection we will prove Car¬ 
leman estimate for a strongly damped hyperbolic equation with variable coefficients. We consider the fol¬ 
lowing hyperbolic equation 

(2.14) (2^ + A)Ar;-div(A*(x)V(9tx) = / in Q = 12x(0,T). 

Many physical phenomena are properly described by the strongly damped hyperbolic equation as (12.14b 
such as equation of this type can be considered as a class of linear evolution equations governing the motion 
of a viscoelastic solid (for example, a bar if the space is M and a plate if the space is ) composed of the 
material of the rate type. 

The following Theorem is a weighted Carleman estimate with a second larger parameter for the strongly 
damped wave equation (12.14b . 














BIOT CONSILIDATION MODEL IN PORO-ELASTICITY 


7 


Theorem 2.3. There exist two constants 7 * > 0 and C > 0 such that for any 7 > 7 ^., there exists s* = 
^*( 7 ) > 0 such that the following estimate holds 

(2.15) cj \Vvf + \vf + a\Vvtf + J |/p 

+ J \^v{x,to)\‘^ + 7~^ \vt{x,to)\‘^ + \v{x,to)f + cP' |V7;(x,io)|^) e^’^'^dxdt 

for any solution v G H^{Q) to problem f l2.i4D compactly supported in Q and any s > s*. 

In order to prove Theorem 12.31 we need the following two Lemmas: 

Lemma 2.4. For any k G K+, the following estimate holds 

f f w{x,T)dT e^^‘^dxdt<C f \'w{x,t)\‘^ e^^'^dxdt, 

JQ Jta Jo 


'Q 

for any w G L?‘{Q). 

Proof By the Cauchy schwartz inequalities, we obtain 

rt 2 


/ ere 

Jq 

Using the fact that, 
we get. 


2S(p 


[ w{x,T)dT dxdt < f a{t — ( f \w{x,t)\‘^ dr) dxdt 

Jtf) Jq ^ Jto ^ 






f a{t — ( f \w{x,t)\‘^ dr] dxdt =—— f [ \w{x,t)\‘^ dT\ dxdt 

Jq ^ Jto ' W Jq Ot V ) 

^ J \w{x,t)\‘^ dr^ _|_ g2s(/j(a;,o) ^ J |y;(2;, r)|^ ^ dx 

+ — /" e^^'^\w(x,t)\'^ dxdt 


1 

'Ip 


Psip{x,T) I 


< 


— f e‘^^'^\w{x,t)\‘^ dxdt. 

4p Jq 


Then 

(2.16) 


f f w{x,T)dT dxdt < C f \'w{x,t)\'^ dxdt 

Jq Jto Jq 


Let A: > 0, and let w{x, t) = a^{x, t)w{x, t), by (12.161) . we obtain 


fie 


,2s(p 


>Q 


/ a^{x,T)w{x,T)di 
Jto 


dxdt < C f 

cr^(x, t)w{x, t) 

Jq 



dxdt. 


Moreover, using the fact that cr(x, f) < cr(x, r) if fo < t < t or f < r < to. we deduce that 


(2.17) [ a^'^+\x,t)e 

Jq 

The proof is complete. 


,2s(p 


w{x,T)dT 


'to 


dxdt < C f a^’^{x,t)e‘^^'^\w{x,t)\'^ dxdt. 

Jq 


□ 
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Lemma 2.5. Let 6 > 0, T > 0. There exists a constant C = (1 + 5‘^T‘^e^^'^) > 0, such that the following 
estimate holds 


dt<C 


tto 


h{T)dT 


Jo J to 

for any h G L^([0, T]). 

Proof For t > to, denote 

fit) = [ e^'^h{T)dT, and g{t) = e~^^f{f). 
dto 

Using the fact that, g'{t) = h{t) — dg{t), we deduce that 

g{t) = / g'{s)ds = / h{s)ds — d g{s)ds. 
J to J to J to 

Then, 

Is-WI < 

Applying Gronwall’s Lemma, we get 


dt, 


f h{s)ds + d f |( 7 (s)| ds 
J to J to 


l5'(i)l < 


h{s)ds 


fto 


+ 5 


fto 


h{s)ds 


Uq 




Then, 


[ \git)\‘^dt< f f h{s)d. 
Jo J 0 J tn 


0 -'to 


dt + d^e^^^ 


rT , rT 


< 


[ [ h{s)d. 

Jo J to 


JO ^ J to 

dt + d'^Te^^^ 


h{s)ds 


to 

*2 ,.T 


to 


dr I dt 


dr 


tto 


h{s)ds 

< (1 + [ [ h{s)di 

Jo J to 


dt. 


For t < to, denote 


fiO = Jt J'"h{'r)dT, and g{t) = e f{t). 


Using the fact that, g'(t) = —h{t) — dg{t), we deduce that 


fto fto fto 

g{t) = — / g'{s)ds= / h{s)ds + d / gis)ds 

Jt Jt Jt 


From the same argument, we conclude 

rT 


f |p(f)|^ df < (1 + f f h{s)d. 

Jo J 0 J to 


dt. 


The proof is complete. 

Now, we give the proof of Theorem l2.3l 
We denote by Pi the parabolic operator given by: 

(2.18) ^i(^) = ut — div(A*(x)Vti). 


□ 
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Then (12.141 ) can be written as 

(2.19) ’P(^) := d'i'^ — (2/U + A)A?; — div(A*(x)Vc)ti;) = Pi{vt) — (2^ + A)Ai; = / 
We deduce that 

(2.20) Pi(ut) = / + (2/r + A)Au. 

Applying the parabolic Carleman estimate give, by Lemmato vt, we obtain 


(2.21) C'y f (a\Vvt{x,t)\‘^ + a^\vt\'^']e^^^dxdt< f \f\^e‘^^'^dxdt-\-C f \Av\^ e^‘^‘^dx dt 
Jq^ ' Jq Jq 

Applying Lemma |A4| to vt, with k = 3/2, we deduce 


/ / Vt{x,T)dl 

J Q to 


e^^'^dxdt < C 


[ a^\vt{x,t)\ 

JQ 


^ e^^'^dxdt. 


Then, we have 
( 2 . 22 ) 


f |u|^ < C f dxdt + C f \v{x,to)\‘^ e^^'^dxdt. 

J Q J Q J Q 


From the same argument, we conclude that 


(2.23) 


[ a^\Vv\'^e^^‘^dxdt<C [ a e^^^^dxdt + C [ \Vvix,to)\‘^ e^^^dxdt. 
J Q J Q J Q 


Collecting (12.221) and (12.231) . we get 


(2.24) 


'Q 


{a^ Iv]"^ + a'^\Vvf)e'^^‘^dxdt < C f (u^ + fi | 

Jq 

+ C f \v{x,tQ)\'^ + \Vv{x,to)\'^)e^^^dxdt. 

Jq 

Inserting (12.211) in (12.241) . we obtain 

(2.25) f \vf + a'^ \Vv\‘^)e'^^'^dxdt < f \ffe'^^‘^dxdt 

Jq Jq 

+ C'^~^ f I Au|^ dt + C f \v{x,tQ)\'^ + \Vv{x,tQ)\‘^)dxdt. 

Jq Jq 

Collecting (12.251) and (12.211) . we get 

(2.26) J |x|^ + |Vxp + cr^ |xi|^ + (T ^dxdt < y {/{^e^^^dxdt 

+ f |Axp e^^'^dxdi + C f \Vv{x,to)\'^ + a'^\v{x,to)\‘^)dxdt. 

Jq Jq 

In order to estimate the second term in (12.261) . denote z{x, t) = Av{x, t), by (12.191 ) z satisfy the following 
differential equation 

zt + S{x)z = - VA*(x).Vxt - /) := G{x,t), 

A*[X) 
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with 5(x) = {2fi{x) + A(x))/A*(x). 

Then using Duhamel Formula, 2 ; takes the form 


(2.27) 


zix,t) = 


'^0 


Let now compute the last integral in (I2.27I) . 
We have 


(^■^{vt{x,T) - VX*{x).Vv{x,t)) - /j 


By integration by parts, we deduce that 
(2.28) 


G(x,r)e‘^(^)(^ = -j^^{{vt{x,t)-V\*{x)Vv{x,t))-{vt{x,to)-V\*{x)Vv{x,tQ))^i 

[ (xt(x,T) - VA*(x).Vr;(x,r))e^(^)(^-*)dr-/ /(x, 

Jto ^ (^) Jto 


,5{x){to-t) 


X*{x) 

Using Lemma [231 we deduce that 


(2.29) J I Ax|^ e^^'^dxdf < C*^ y Ax(x, to)|^ + |u(3^) fo)|^ + |Vx(x, e^^^dxdt 


'Q 


Vt{x,T)dT 


e^^'^dxdt 


e^^'^dxdt + 


'Q 


f{x,T)dT 


z^^'^dxdt^. 


+ J {\vt\‘^ + \^vf^e‘^^^dxdt + 

+ f f Vv{x,T)di 
J Q Jto 

Consequently, we find 

(2.30) J |Vx|^ + (T^ |xt|^ + (T |Vx4p ^dxdf < y \f\‘^e‘^^‘^dxdt 

+ C [ e^®‘^(cr^ |Vx(x,fo)P + \vtix,to)\'^ + 7 “^ \Av{x,to)\‘^)dxdt. 

Jq 

This completes the proof of Theorem [23] 


2.3. Carleman estimate for hyperbolic equation. We recall the following Carleman estimate for a scalar 
hyperbolic equation. As for the proof, we can refer to Bellassoued and Yamamoto and Isakov and Kim 
Ii2^ l25l which gives a direct proof by integration by parts. 

Lemma 2.6. Let /r > 0 satisfy | |V (log /r)| |x — xo| < 1. There exist constants C > 0 and 7 * > 0 such 
that for any 7 > 7 *, there exist s* = ^*( 7 ) such that for all s > the following Carleman estimate holds 

(2.31) cj (T^^lVyp + + 0 -^ < J \ytt - lx{x)Ayf e'^'^'^dx dt, 

for any y E H^{Q) with compact support in Q. 
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2.4. Completion of Carleman estimate for Biot’s system. Let (v, y) be a solution of Biot’s model in 
poro-elasticity (12.101) satisfying (12.1 II) . 

Let v{x, t) = div v(x, t), we apply div to the first equation in (12.101) . we can derive the following equation 
(2.32) utt — (2|U + A)Au — div(A*Vut) — VA*.Vu 4 + piAy + V^i • Vy = divf in Q. 
Applying Theorem l2.3l to (12.321 ). 


(2.33) cj + a^\vf + a\Vvt\‘^+ a^\vt\^y^^^dxdt<-f-^ j |div fp 

+ j |Au(x,to)P + 7“^ \vt{x,to)\‘^ + \v{x,to)f + |Vu(x,fo)P) e^'^'^dxdt 


+ 7 M (|Ay|^ + 

JQ 

Now, let w{x,t) = rot v(x,f), similarly we apply rot to the first equation in (12.101) . we can derive the 
following equation 

(2.34) wtt — yAw = rot f in Q. 

Applying Lemma |2^ of the hyperbolic Carleman estimate, we obtain 

(2.35) C J (T^ |Vt(;|^ + + (T^ < J |rot f|^ df. 

We recall that, the first equation in (12.101) . can be written 

(2.36) vuix,t) - y{x)Av = fo{x,t) in Q, 
where, 

fo{x,t) = i{x,t) + (y + A)Vu + uVA + (Vv + (Vv)'^)V/U + V(A*(x)ut) - pi(x)Vy. 
Applying again Lemmaof the hyperbolic Carleman estimate, we obtain 

(2.37) C J I Vv|^ + |vt(x, f)|^ + (T^ |v|^ ^ e^^‘^dx df < J |f e^®‘^dx df 

+ ^ + \^t\^ + + 1 ^ 1 ^ + |Vy| 2 )e 2 '*^dxdL 

Collecting (12.331 ). (12.351 ) and (12.371 ). we get 


(2.38) 

c J (^<T(|Va;,tv|^ + |Va;,trup + (T |Vu|^ + (T^ |ut|^ + |Vut|^) +cr^(|v|^ + cr|u|^ + |t(;p)^c^^'^dxdf 
< [ (|f|V |Vf|2)e2*^dxdt+ [ ( 7 '^ |Ay| V |Vyp)e^'*^dxdf 

jQ Jq 


+ j ^ |Ax(x,fo)|^+ 7 ^ |xt(x,fo)|^ + |Vx(x,fo)|^) e^'^'^dxdf. 

We recall that, the second equation in (12.101 ) is given by: 


yt{x,t) - Ay{x,t) + y 2 ( 2 :)divvt(x,f) = h{x,t) 


in Q. 
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Furthermore, applying Lemmaof the parabolic Carleman estimate, we obtain 
(2.39) C 7 / (^\Ay\^+ a^\Vy\^+ a^\y\^y^^^dxdt 

J icJ 


< 


f a \h\^ dx dt + f a |div vt|^ e^^'^dxdt 
Jq Jq 


Then we have, 

(2.40) C'y~^ J ^ I Ayp + |Vy|^ + (T^ |?/P ^ (if 


<7 ^ / a\h\'^ dx dt + ^ ^ f a\vt\^e^^'^dxdt. 
Jq Jq 

Inserting (12.401) into (12.381) . we find 
(2.41) 

C J y I Va;,tvp + a \ Vx,tw\'^ + cj^ I Vnl^ + (7^ l^tl^ + (7 I Vvil^ + (7^ |v|^ + cj^ Itup + cT^ |r;|^ ^e^®‘^cix dt 

< [ {\if+ \V{f+-f-^a\h\‘^)e^^‘^dxdt 


'Q 


+ ^\Av{x,to)\‘^ + 'y ^\vt{x,to)f + (T‘^\v{x,to)\‘^+ a‘^\Vv{x,to)y e^'^'^dxdt. 

J Q 

We deduce that 

(2.42) C J ^(7 |Va;,tvp + (7^ |vp + (7^ |r;|^ + (7^ + (7^ iVnp + (7 iVuil^ ^e^®‘^cix(ii 

< [ {\if+ \Vi\'^+-f-^a\hf )e^^^dxdt 

Jq 

+ J ( 7 ”^ |An(x,io)P + 7”^ |^^((2:,io)|^ + |^^(a^,io)|^ + <7^ |Vn(x,io)|^) e^'^'^dxdt. 

Additionally (12.391 ) and (12.421 ). we find that 

(2.43) cj (^C7 |V3;,iVp + (7^ |vp + (7^ |r;|^ + (7^ + (7^ iVnp + (7 iVxil^ 

+ |A(/|^ + (7^ |Vyp + |y|^ ^e^^'^cfxcii 

< [ (|f|V |Vf|V7"V|/i|2)e2"‘^cix(ii 
Jq 

+ J (^7~^ |Ax(x,io)|^ + 7”^ |^^((a;,io)|^ + <7^ + <7^ |Vx(x,io)P) e^'^'^dxdt. 

The proof of Theorem 12. ll is complete 

3. Proof of the main results 

This section is devoted to the proof of Theorem 11.11 and Theorem 11.31 The idea of the proof is based on 
the Carleman estimate method. 
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A usual methodology by Bellassoued and Yamamoto |4], yields an estimate of a source term A*(a:) and the 
two spatially varying density. 

3.1. Preliminaries estimate. 

Lemma 3.1. Let oj be an open subdomain of LI with regular boundary dco D F. There exists constants 7*, 
s* and C > 0 such that for any s > and any 7 > 7 * the following estimate holds: 


(3.1) 


/ 

J U) 


a 


^ lup e^^'^dxdt < C 


u}x{0,T) 


lujx(0,T) 


cj^ |Vup e^^'^dxdt 


for any v G H^{uj x (0, T)) such that v{x, t) = 0 on du x (0, T). 

Proof We multiply Vu by and using the divergence theorem, we obtain 

f VV.{yp)v dx = — f vdiv{{V(f)v e‘^^‘^)dx 

J UJ J ijJ 

— ~ [ \v\‘^ Aif e^^‘^dx — 2s f |u|^ | 

J LJ Juj 

— f V V.yip) V dx . 

J ijJ 

Therefore, 

2 f aVv.yd)v e‘^-^‘^dx = -2 f \vfyb\^ e'^^^dx 

Juj J LJ 

— f cr |u|^ Ad — 7 f u |r;|^ |Vdp e^^'^dx. 

J UJ J (jJ 

Taking 7 > 7 * and s > s* sufficiently large, we obtain for any e > 0 

(3.2) C f a^\v\‘^ e^^'^dx <Cs f iVxpe^'^dx + i 

J (jj J (jj 

Integrating in (0, T) and taking e small we obtain 


a^lvl^e^^^dx. 


(3.3) 


C 


a'^\v?e^^^dx < a 


Jujx{0,T) 

Applying the last inequality to a'^v we obtain 
(3.4) [ a^\v\‘^e^^^dx<C 

Jujx{0,T) 

for any 7 > 7 * and s > s*. This completes the proof. 

Hence, by Lemma [TTI we obtain the following Lemma. 
Lemma 3.2. Let (v,y) G H‘^{Q) x satisfying 


[ yvfe'^^'^dx. 
Ju}x{0,T} 


[ a^yvfe'^^'^dx, 

Ju}x{0,T) 




Vti 

- A^ 7 V - V(A*(x)divvt) + 

gi{x) Vy 

II 

G 

Q, 

(3.5) 

< 

Vt - 

Ay + Q 2 {x) divvt = h 


(x,f) 

G 

Q, 



V = 

< 

II 

0 


(x,f) 

G 

s, 

and 








(3.6) 

d{sr{x, 0 ) 

= 

v{x,T) = 0, y{x,0)=y{x 

T) = 0 

for all X 

G 

n, 


□ 
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There exist positive constants 7 * and C > 0 such that, for any 7 > 7 * we can find s^, and D, the following 
inequality holds 


(3.7) J |V 3 ;^fVp + |v|^ + |divvp + |divvip + cr^ |Vdivvp + a |Vdivvt|^ 

+ |Ay|^ + |Vyp + < C J ^7“^(T |/ip + cr^ |f|^ + |Vf|^ 

+ C'e^"(||v||^4(^x(o,T) + l|divv(-,to)|lH2(n) + l|divvi(-,to)|li2(Q)) 


for any s > s*. 

Proof Let u' C uj such that du' D L. In order to apply Carleman estimate, we introduce a cut-off function 
X satisfying 0 < x < 1> X ^ x = 1 in ^\uj' and suppx C Tl. 

Put 


v(x, t) = x{x)^{x, t), y{x, t) = x{x)y{x, t), 
Noting that (v,y) £ H'^iQ) x satisfying 


(3.8) 


with 


where 


I vtt - -V(\*divvt) + giVy=? (x,t) € Q, 

[ yt- Ay + g 2 divvt = g (x, t) E Q, 

Supp(v(-,f)) C n Supp(y(-, f)) C fl, VfE(0,r), 

f = x{x)f{x,t) - [A^,a,x]v - V(A*Vx.vt) - A*(Vx)divvt -f piyVx, 


g = x{x)h{x, t) - 2VxVy - yAx + Q2^X ' v* 

and 

[A;,,A, x]v = A^,a(xv) - xA^,xsr. 

Noting that (v, y) satisfies (13.81) . then we can apply the Carleman estimate for Biot’s system (12.121) to (v, y), 
we obtain 


(3.9) 

C 


'{n\ui')x{0,T) 

+ 


I Va;^tv|^ -f |vp -|- cr^ |divv|^ -|- |divvt|^ -|- |Vdivv|^ -|- a | Vdivv^l 
\Ay\^ + a‘^\Vy\^ + a'^\y\^'^dxdt< f a \gf + 

J Q 


dxdt 


J ^ IAdivv(x, fo)|^ + 7 ^ |divvt(x, fo)|^ + |divv(x, fo)|^ + IVdivv(x,fo)|^^ 
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By a simple computation, we get 


(3.10) I e2^^(7-V|p|2 + 


Vf 


)dxdt <cj^ ( 7 -V \xh\^ + \xif + \V{xi)\^)dxdt 
+ C ( |(Ax)v |2 + |Vx|' |Vv |2 + |V[A^,A,x]v| 2 )e 2 ^‘^dxdt 
+ ( |V(A*Vx.vt)|2 + |A*(Vx)divvi|2 + |V(V(A*Vx.Vi))|' + |V(A*(Vx)divvi)|2 )e2*^dxdt 

+ C [ (\iAx)yf + \Vx\^\Vy\^ + \Vx\"\y\^)e^^^dxdt. 


'Q 


Since SuppAx, SuppVx C tu' we obtain 

(3.11) 


'{OV)x(0,T) 


|v|2 + |divv|2 + |divvtp + iVdivvp + a iVdivvtp 

+ \Ay\^ + |Vy |2 + |y |2 )dxdt < C (x'V \h\^ + jfp + |Vf |2 )dxdt + ||v|| 1 , 4 (,^(o,t)) 

+ c( [ a‘^ \y\^ dxdt + f a'^\Vyf e^^'^dxdt) 

^ Ox(0,T) / 


a;'x(0,T) 


+ (x ^ IAdivv(x, to)|^ + 7 ^ |divvt(x, to)|^ + |divv(x, to)|^ + IVdivv(x, e^^'^dxdt. 


IQ 


Let xi be a cut-off function satisfying 0 < X1 < 1. Xi G Xi = 1 in cu' and Supp(xi) C oj. Let us 

consider z{x,t) = Xiix)y{x,t) & H^{uj x (0,T)) and z{x,t) = 0 for all (x,f) & duj x (0,T), so that by 
Lemma l3dl we have 

(3.12) f \y\'^ ae^^'i^dxdt < f \z\'^ e^^^^dxdt 

Joj'x(0,T) Ju}x(0,T) 

< 


C [ a‘^\Vyfe^^^dxdt + c[ \y\^ e^^^dxdt. 

Ju}x{0,T) Jq 


Furthermore by the first equation of (13.51 ). we have 


(3.13) f (T^ |V?/|^ e^^'^dxdf < C'e'^* ||v||^3. .pN-j-h f |fp e^^'^dxdf. 

Juj'x{0,T) ^ Jnx{o,T) 

Inserting (13.131) in (13.111) . we obtain (13.71) . This completes the proof of the lemma. 


□ 


Henceforth we fix 7 > 0 sufficienfly large. By Ng^^p we denofe fhe quantify 

(3.14) A,,^(v, y) = |Vx,iv|2 + |v|2 + |divv|2 + 

+ |Vdivv|^ + s iVdivv^l^ + \Ayf + |Vyp + |yp 'jdxdt. 

We infroduce a cul-off function rj satisfying 0 < rj < 1, y £ C°°(M), r/ = 1 in (2e, T — 2e) and 
Supp( 7 ) C {s,T — s). Finally we denofe 

v = r/v, y = vy- 
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Setting do = < 5 / 2 ) 7 ^ we have 

max?/)(x, t) < do) i G (0,2e) U (T — 2e, T). 

x&Q, 

Lemma 3.3. There exist three positive constant s^, C > 0 and D such that the following inequality holds: 
C'A^s,^(v,^< j (s|/ip + s2|f|2 + |Vfp)e2*^dxdt 

+ (l|v|lH4(^x(o,r)) + l|divv(-,io)|lH2(n) + l|divvt(-,to)|li2(Q)) 

+ (||v|Ih1{o,T;H1(o)) + \\y\\l^(Q)) 

for any s > s* and any (v, y) G H‘^{Q) x satisfying 

wtt - A^,av - V (A*divvt) + giVy = f (x, f) G Q, 
yt - Ay + £>2 divvt =/i {x,t)£Q, 

V = 0, y = 0 (x, t) G S. 

Proof Let (v, y) G H'^{Q) x put 

v(x, t) = y(t)v(x, i), y(x, t) = y{t)y{x, t). 

Noting that (v,y) G H^{Q) x satisfies 


vti - A^_av - V(A*divvt) + yi Vy = yf + rjtt^ + (x, t) G Q, 

(3.15) yt - Ay + y 2 divvt = y/i + yt(y - ydivv) {x,t)eQ, 

V = 0, y = 0 (x, f) G S, 

and applying Carleman estimate (13.71 ) to (v, jf), we obtain 


CNs,^(y,y)< j (s|d|2 + s2|f|V|Vf|^)e2*^dxdf 

J Q 

+ (l|v|lH4(^x(o,r)) + l|divv(-,fo)|lH2(n) + l|divvt(-,fo)|li2(Q)) 

+ f (s"(|y«P + |y4p)(|yp + |v|2 + |v4|2 + |Vv|2 + |Vvt|2)e2'*^dxdf 

JQ 


for any 7 > 7* and s > s*. Since Supp(ytt), Supp(yt) C (0, 2e) U (T — 2e, T), we obtain from (12.301) 
f {s^{\vtt\^ + \m\^){\y\^ + |vp + Iv^p + |Vv|2 + |Vv4|2)e2'*^dxdf 


'Q 


< Cq^p^'^OS ( ||.^||2 

s OS e \^l|V||j^l(0,T;Etl(O)) 

This completes the proof of the lemma. 

Lemma 3.4. There exists a positive constant C > 0 such that the following estimate 
/ \z{x,tQ)\^ dx < C / (u l 2 ;(x,f)|^ + | 2 ;t(x,f)|^)dxdf 

Jcl j q 

for any z G 77^(0, T; L^(ff)). 


+ 


□ 
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Proof. By direct computations, we have 

rto ^ 


j ri^{tQ)\z{x,tQ)\^ dx = J ^(y 'n'^{t)\z{x,t)f dx^dt 

[ [ 'n'^{t)z{x,t)zt{x,t)dxdt + 

Jo Jn 


rto 


= 2 


10 Jn 


r]t(t)rj{t) \z{x,t)\ dxdt. 


Then, we have 


/ \z{x,to)\‘^ dx < C / {a\z{x,f)\^ + a ^\zt{x,t)\^)dxdt. 

J Q, J Q 

This complete the proof of the lemma. □ 

Finally, let <p{x,t) be the weight function defined by 

= :=p{x)a{t), 

where, p{x) and a{t) are defined by 

(3.16) p{x) = = d,\/x G 0 and a{t) = < l,Vf G (0,T). 

3.2. Proof of the stability in determing A*. We prove now Theorem ll.il To this end we use the global 
Carleman estimate (12.121) . 

Consider now the following system 

utt - A^^au - V(A*(x)divut) + qi{x)V9 = 0 {x,t) G Q, 

(3-17) 6t - A6 + g 2 {x)divut = 0 {x,t)GQ, 

u = 0, 0 = 0 (x, t) G S. 

Henceforth, for simplicity, consider the following functions: 

u = u(A*, pi, P 2 ), u = u(A*, pi, P 2 ), 9 = d{y,Qi,Q 2 ), 9 = 9(x*,qi,Q2) 

and 

v = u —u, y = 9 — 9, f = X* — X*. 

Then by (13.171 ). we easily see that 

vtt - A;,,av - V(A^divvi) + piVy = V(/(x)divut) (x, t) gQ, 

(3.18) < yt-Ay + Q 2 divvt = 0 (x, t) G Q, 

V = 0, y = 0 (x, t) G S. 

In this subsection we discuss a linearized inverse problem of determining A*. We assume that the assump¬ 
tions (A.l) and {A.2) holds true, then our inverse problem is identification of /(x). 

Let 

v^ = dlv, y' = dly, i = 1,2,3. 

By a simple calculation, we obtain for t = 1, 2,3, (v*, y*) satisfies 

vlt - A/,,av* - V(A*(x)divvJ) + pi(x)Vy* = V(/(x)divuJ) (x,t) G Q, 

(3.19) yl - Ay^ + g 2 (x)divvl = 0 (x,t) G Q, 

V® = 0, y* = 0 (x,t) G S. 
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Let V = divv, we apply div to the first equation of the system (13.181 ). we get 

(3.20) A(/(x)divui)(x,fo) = 

[vtt - (2// + A)Au - /\{X*{x)vt) + Qi{x)^y + Vpi(x).Vy)(x, to) = 

We have f{x) = 0 and V/(x) = 0 on F and 

|divut(x,fo)| > e, 

then by the elliptic Carleman estimate, we deduce that 

J (^s^\f{x)f+ s\Vf{x)f+ s-^\Af{x)f^ < C J \k{x,to)f 

<C [ \vtt{x,to)f 

Jo, 

+ (l|divv(-,to)||^2(t2) + ||divvi(-,to)||^2(o) + ll2/(-,fo)|lH2{n)) ■ 

Applying lemma |A4] to z{x,t) = vuix, t) = we get 

J \vtt{x,to)\'^ < C J (^s\vttix,t)f + s~^\vttt{x,t)f'^ 

<C j {s |divvi(x,f)|^ + |divv2(x,f)|^) e^^^^^'^'^dxdt 

We deduce that 

(3.21) J (^s5|/(x)|^ + s3|V/(x)P + s|A/(x)|2) 

(^||divv(-,to)||^ 2 (o) + ||divvt(-,fo)||^ 2 (o) + 

Now, applying lemma [F3] to (v®, y*), for f = 1,2, we obtain 

(3.22) CNs^^{^\f) < (^s‘^ |V(/(x)divUt)|^ + | A(/(x)divUj)|^^ e^^'^dxdf 

+ (lK1li4(^x(o,T)) + l|divv'(-,fo)t2(f,) + ||divv®(-,fo)||'2(^)^ 

+ e ° (^||v ||j^i(o,T;/ti{n)) + \\y IIl 2 {q)^ 

Inserting (13.221 ) into (13.211 ). and using Assumption (A.2), we obtain 

J (s® \f{x)f + |V/(x)|2 + s |A/(x)12) 

Y1 (ll'^1li4(^x(0,T)) + I|divv'(-,fo)t2(o) + ||divv®(-,fo)t2(n) + ll2/(-,io)|ll^2(n))+CsV'='°*Mo 


< e 


i=l,2 
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We deduce that 

s® |/(x)|2 + |V/(x)p + s |A/(x)|2) 

< (^I|v|Ih6(^x(o,t)) + lly(->*o) 11^2(0) + J^||divaiV(-,to)||5^2 (o)^ + CsV'^°"Mo 

< (l|v|lHa(,,(o,T)) +A'to(v,2/)) + Cs 

Finally, minimizing the right hand side with respect to s, we obtain: there exist k E (0,1) such that 

(3-23) ll/lliT2(o) — C* (^||v||^6(^^x(o,t) +y)) • 

The proof of Theorem [TTT] is completed. 

3.3. Proof of the stability estimate in determining gi and Q 2 . For simplicity, we set 

U = u(A*,Pi,P2), U* = u(A*,Pi,P2) 

and 

e = eix*,Qi,Q2), 9* = 9{x*,Qi,h). 

Let (u, 9) satisfies the following equation 



(3.24) 

and (u*, 


/ 

utt{x,t) - A^^xu(x,t) - V(X*divut(x,t)) + giV9(x,t) 
< 9t(x,t) — A9(x,t) + g 2 divut(x,t) = 0 
u(x,t) = 0, 9{x,t) = 0 

9*) satisfies fhe following equafion 


0 in Q 

in Q 
on S 


(3.25) 


Lef 


/ 

u*^(x,t) — A^^;,u*(x,t) — V(A*divu*(x,t)) + piVd*(x,t) = 0 in Q 
< 0*(x,t) — A0*(x,t) + p 2 divu*(x,t) = 0 in Q 

u*(x,t) = 0, 9*{x,t)=0 on S 

v = u — u*, y = 9 — 9* 


p{x) = gi{x) - gi{x), and q{x) = g 2 {x) - g 2 {x), 

where (u, 9) satisfies (13.241) and (u*, 0*) satisfies (13.251) . Then, by a simple calculation, we have 


(3.26) 


wtt - A^,av - V(A*divvt) + giVy = p{x)V9* {x,t) E Q, 

yt- Ay + P 2 divvt = q(x)divu5‘ (x, t) E Q, 

V = 0, y = 0 (x, t) E S. 


In (his subsecfion we discuss a linearized inverse problem of defermining gi and p 2 - We assume thaf (he 
assumptions (A.l) and (^4.2) holds (rue, (hen our inverse problem is identification of p{x) and q{x). 
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For i = 1, 2, 3, we denote v* = v, y® = dly, where (v, y) satisfies (I3.26I ). Then, we have 


(3.27) 


- A^,av' - V(A*divv®) + yiVy' = p{x)Vdie* {x, t) € Q, 


yi - Ay® + y 2 divvj = y(x)div5iut* 

V® = 0 , y* = 0 


(x,f) E Q, 
(x, f) E S. 


We apply lemmato (v®, y®), for z = 1, 2,3, we obtain 

(3.28) J (^s\q{x)divdiu*tf + s^\p{x)Vdte*f + \V{p{x)Vdi9*)f^ e'^^'^dxdt 

+ (lK1li4(^x(o,T)) + l|divv'(-,fo)t2(f,) + ||divV®(-,fo)||'2(^)) 

+ Cs e ° (^||v ||j^i(o,T;/ti{n)) + lb IIl2(q) 

where 

v® = y(f)v\ ^ = v{t)y\ i = 1,2,3. 

In terms of Assumption (A.2), we obtain 

(3.29) C'As^y,(v®,^) < / |g(x)|^ + |p(x)|^ + |Vp(x)p^ 

J Q 


+ e 


Ds 


+ ||divv®(-,to)||^2(t2) + l|divvK-,to)||i2(t^)) + C's^c^'^o^Mq. 


ll/t4(^x(o,r)) 

In the following we give two lemmas which will be used in the proof of the stability theorem of determin¬ 
ing of the two spatially varying coefficients. The first one is a Carleman estimate for the first-order partial 
differential equation: 

We consider a first order partial differential equation: 

3 

P{x,D)i = E 7j(x)(9jf-f 7o(x)f, xEfi, 
i=i 


(3.30) 


where 

(3.31) 
and 

(3.32) 


7o(x) E W ’°°(f2), 7(x) = (7i(x),72(x),73(x)) E (W 


| 7 (x).(x — xo)| > Co > 0, on 12, 
with a constant cq > 0. Then we have the following Lemma 

Lemma 3.5. In addition to d3.3iD and d3.321) . Then, there exist constants s* > 0 and C > 0 such that 

(3.33) |fpdx < C [ |P(x, Z2)f|^ dxdf 

Jq Jq 

and 

(3.34) 


' [ iVfl^dx <C [ e2"‘^(*'*)(|P(x,I2)f|2 + \V{P{x, D)f)\‘^)dxdt, 

Jq j ft 


for all f E satisfying f(x) = 0, Vf(x) = 0, x E T and all s > s*. 
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Proof. Let denote = l{x) ■ (x — xq). We multiply the both sides of (13.301 ) by f (x)C(x)e^*^(*) and 
using the divergence theorem, we obtain 

(3.35) f P(x, Z))f(x)f(x)C(x)e^*^*'*^(ix 

Jn 

= f Vf(x). 7 (x)f(x)C(x)e^^^^^^(ix + f ^Q{x)\i{x)\^ dx 

Jq J q. 

= — f f(x)div(7(x)f(x)C(x)e^*^*'^^)(ix + f '^o{x)f{x)\^C{x)e^^'^^^'^dx 

Jq, J d. 

= — f |f(x)|^ div(C(x) 7 (x))e^^^*^*^(ix — f 7 (x) • Vf(x)C(x)f (x)e^^^*'*^dx 

Jq J q 

— 2sf |f(x)p C^(x)e^*^*'*^(ix + f 7 o(x) |f(x)|^ 

j q. 

Using (13.321) . and the fact that 

7 (x) • Vf(x) = P(x, D)i - 7 o(x)f(x) 

We deduce, in terms of (13.351 ) and the Cauchy-Schwartz inequality, 

2 c^s [ |f(x)|^e2'*^(^)dx < 2s f \i{x)f 

Jq, j q. 

< 2 [ |P(x,L>)f| |f(x)|e2^^(^)dx+ [ | 7 o(x)C(x)| |f(x)|2e2*^(*)dx 

Jq Jq 

< c [ |P(x,D)f| |f(x)|e2*^(^)dx + U [ |f(x)|2e2"^(^)dx. 

Jq, j d. 

Then for large s, we get 

s[ |f(x)|2e2^^(^)dx < C [ |P(x,U)f| |f(x)|e2"^(^)dx. 

J D. J Q. 

On the other hand, for all small e > 0 there exist a constant > 0 such that 

(3.36) [ |P(x,P)ff(x)C(x)|e2*^(^)dx < —[ |P(x, P)f ^ + es / |f (x)|2 

Jo. ^ Jo. Jq, 

In terms of (13.361) . we have 

s [ |f(x)pe^'*^(^)dx < 

Jn 

Moreover, for all small e > 0, we get 


[ \P{x,D)i\ 
Jn 


2^2spix)dx^SS [ \i{x)fe^^P^Pdx. 
Jn 


: [ |f(x)| 
Jn 


^e^^PPdx < — [ 
s Jn 


Consequently, 


Since 


[ |f(x)l 

Jn 


2 ^2sp(x)^^ < (J 


\P{x,D)i\^ e^^PPdx. 

[ \P{x,D)ife‘^^PPdx. 

Jn 


P(x, D)dj{ = djPix, P)f - idjP{x, P))f, dji |r= 0, 
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we apply (13.331 ) to dji, for j = 1,2,3, we get 

^2 f \djif dx < C [ e‘^^'P^^’^^\P{x,D)djif dxdt 

Jn Jn 

<C [ \djP{x, D)f - {djP{x, D))f\‘^ dxdt 

Jn 

<C [ + \djif)dxdt + C [ \djP{x, D){\‘^ dxdt 

Jq Jq 

<C [ e^^^^^’*\\P{x,D)i\‘^+ \djP{x,D)i\‘^)dxdt + C [ \Vif dxdt. 

Jo. Jo. 

Therefore, 

g2 f e 2 «p(^) iVfpdx < C [ e‘^^^^^'^\\P{x,D){\^+ \VP{x,D){\‘^)dxdt 

Jo, Jn 

+ C [ |Vf|^dxdt. 

Jn 

For sufficiently large s, we can complete the proof of the lemma 1331 

□ 


The second one, is 
Lemma 3.6. Let (v*, y®) satisfy 


(3.37) 


vJt - A^_av® - V(A*divvJ) + piVy® = f® (x, t) G Q, 
< _ Ay® + P2divvj = y® (x, t) G Q, 

V® = 0, y* = 0 (x,f) G S. 


Then there exist constants C > 0, s* > 0,/or i = 1, 2, 3, such that 


(3.38) 


s [ e^^^^^’^°^\fix,to)fdx< 

Jn 


>Q 


^2sip dxdt + CN„{v\y^) 


+ se 


2sdQ 





(3.39) 

s2 [ e2*^("’‘°)|v®(x,^o)|'dx<c(fV,,(v^^) + s-2^^,(v/,y/)). 

t/ 

(3.40) 

^2 f ^2s^(x,to) |divv®(x,fo)|^da: < C (^fV^(v%y®) + s~'^Nfy^\,y\)^ 

J 

and 


(3.41) 

[ Vdivv®(x,to) ^dx < CNfyvfyi), 

Jn 

for all s > s*. 
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Proof. Applying Lemma l3A] to z^{x,t) = t) = i) and by the second equa¬ 

tion of (13.371) . we obtain 


f \y'^{x,to)\^ dx < s f |^(x, t)|^ dxdf + s ^ f \yf{x,t)\^ dxdt 

J ^ J Q J Q 

<s f |^(x, f)|^ dxdf + f \g^{x,t)\‘^ dxdt + s~^ f | A^(x, f)|^ dxdf 
Jo Jq Jq 


+ S 


Jq 


'Q 

^2s<p 


divvl 


dxdt + e °(||v ||j^l(o,T;iTl( 0 )) 11 ^ IL 2 (Q)) 


Then, 


Cs [ \y\x,to)f dx < [ \g\x,t)f dxdt 

Jn Jq 


+ 


I e2*^(|Ar| 

L 


Q 

2 2 I —)• 12 

+ S^\^\ + 


divv/ 


dxdt + se'^^‘^°{\\v^f 


H^{0,T-,H^{n)) 
2 


+ y 


Il 2 (q) 
|2 


< \9\x,t)\ dxdt + CN^ (v*, y*) + ( 11 v* 11 (q)) + 112 / 1 |l 2 (q))• 

Applying Lemma l34l to z^{x,t) = f) = v*(x, t), we obtain 


dxdt 


g 2 f \yr\x,to)\^ dx < \^\x,t)f dxdt + s [ v/(x,f) 

Jo Jq Jq 

Applying Lemma lL4l to z^{x,t) = g{t)e'^^'^^^J^divv^{x,t) = e^*‘'’*^*’*)divv*(x, t), we obtain 


„2s(/j 


divv/(x,t) 


dxdt 


g 2 f to)|^ d® < f |divv*(x,f)|^ dxdf + s f 

J ^ J Q J Q 

<c(A,(v\^) + s-2Ar,(vAy/ 
Finally, applying Lemma lL4l to z^{x,t) = g{t)e^^‘^^^J^'Vdivv^{x,t) = e^^‘^(*’*)Vdivv*(x, f), we obtain 


/ 

Jo 


^2stp{x,to) 


^2sip 


Vdivvt*(x,f) 


dxdt 


|Vdivv*(x, fo)|^ dx < s f |Vdivv*(x,f)|^ dxdf+s ^ f 

Jq Jq 

This complete the proof of the Lemma. □ 

Now we complete the proof of Theorem 1 1.31 let v = divv. Then we have 
(3.42) div(p(x)V 6 '*(x,to)) = ivu - {2fJ. + X)Av - A{X*vt) + div(piVy)) (x,to) = -^(a^,fo)- 

Note that by the assumption (A.2), we have p{x) = 0 and Vp(x) = 0 on T and 


|Vr(x,fo)| >e. 
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Furthermore the assumption (11.71 ) imply that 

(3.43) [ \F\‘^ < C [ |d ivv^ 

Jq Jn 

+ C [ \Av{x,to)fe‘^^‘^^^’^°'>dx + C [ |A(A*r;t(x,to))|^e2*'’^'^’‘°)dx 
Jq Jq 

+ C [ |div(piVy(x,to))|^e2"^^^’*o)dx. 

Jq 

Using (13.401 ) in the Lemma |3]6j we get 

(3.44) f < C ) ,y^)) 

Jq 

+ Ce^^ (^||divv(.,to)||^ 2 (t 2 ) + ||divvt(.,to)||^ 2 (t 2 ) + \\yi;to)\\H^Q)) ■ 

Similarly, we have 

(3.45) [ \VFfe'^^‘^‘'^’^°'>dx<c[ |Vdivv2(x,to)|^ 

Jq J q 

+ Ce^^ (^||divv(.,to)||^3(Q) + ||divvt(.,to)||^3(Q) + ||2/(-,io)|lH3(n)) • 

Using (13.411 ) in the lemma IX6l we get 

(3.46) 

^ < CJV„(v^5")+Ce'’• (||divv(.. i„)||^.,n, + ||divv,(., + ||!,(., MIlLcn) 

By the second equation of (13.261) . we have 

(3.47) ytix,to) = q(x)divuj*(x,fo) + Ay(x,fo) - P 2 divvt(x,to)in ^ 

Then, 

(3.48) C\q{x)\ < \q{x)\ |divu*(x,to)| < |Ay(x,to)| + |divv^(x,fo)| + ■ 

Moreover, 

(3.49) Cs [ |q'(x)|^dx < s f e^®‘^*'^’*°^(|Ay(x,to)|^ + |divv^(x,fo)|^ + 

Jd. Jq, 

< S [ + s |divvi(x, fo)|'dx + ||2/(., to)|lH3(n) 

Jq Jq 

Then, using Lemma |3^ we get 

(3.50) s [ \q{x)\^ dx < C {N^{v\y^) + 

Jq 

+ + ||y^||i2(Q)) + ll2/(->^o)|lH2(n) ■ 

Then again, by (13.421 ). we observe that 

F{x,to) = Vp(x) • V0*(x,fo) + p{x)A6*{x,tQ), x € fi. 
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So, we need Carleman estimate for the first-order partial differential equation given by Lemma [331 then we 
have 


(3.51) 
and 

(3.52) 


Jq Ju 


f |Vp(x)|^dx < C [ e2"^(^’*)(|Fp + |VFp)dxdf. 

Jq, J d. 

Moreover, using (13.441) . (13.461) and (13.501) 

J g2s¥j(a;,to) ^g4 ^p(^x)f + |Vp(x)|^ -f s |q(x)|^^ dx 

< J (^s2|F(x)|V|VF(x)|Vs|q(x)|^)dx 


-f 


< C (iV.(v\ y^) + + s-^N„{w\f)) 

(^||divv(.,fo)|lH3(n) + l|divvt(.,fo)||^3(n) + l|y(-, io)||^3(n)) +Cse^^'^°M^. 


Furthermore, 


J g 2 s¥>(a;,to) ^g4 |p( 3;)|2 _|_ ^2 ^ dx < C J s\qf e^^'^dxdt 


+ Ce®*E IF 


ltT4(^x(o,r)) 


+ ||divv*(-,to)||?^3(n) + l|divvK-,to)||?^3(n) + 


m{n) 


2 = 1 




for all large s > 0. 

Using the Lebesgue theorem, we get as s —)• oo 

Cs [ |g|2g2s9=(^4o) r ^ g f |^| 2 g 2 sv(x,to)^ 2 .^ 

Jn Jo Jn 

Then, we get 

f ^ 2 s^ix,to) |'g4 ^ ^2 |Vp(x)|2 + S |q(x)p) dx 

u 

< C'e'^® ( l|v||^7(t^x(o,T)) + II2 /(')^o)|Ih3(q) + ^ ||divc)^v(-,fo)||j;^3(Q)^ + CMq , 


2 = 1 


for all s > s*. 

Using (12.71) and (12.81) . we obtain 


^2sd 


In 


^ + |Vp(x)|^ + | 9 (x)|^) dx < C'e^^(||v||^7(^^(o^.p)) + Mt^{w,y)) + 








26 


M. BELLASSOUED AND B. RIAHI 


for all s > s*. 

Moreover, 

(3.53) ^ (|p(x)|2 + \Vp{x f + |(z(x)|2) dx 

for all s > s*. 

Besides, 

(3.54) J (|p(x)|V|Vp(x)p + |Ap(x)|2)dx 

< +Mi„(v,2/)) 

for all s > s*. 

Finally, minimizing the right hand side with respect to s, we obtain: there exist 6 G (0,1) such that 
IIfIIhf^) + lkllL2(n) ^ c* (ll’^ll/t7((^x(o,T)) + ^to(v,y)) • 

The proof of Theorem 1 1.3 1 is completed. 
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